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Recent progress of lattice investigations in thermal physics is summarized in this 
contribution. Hadronic spectral functions can be reconstructed from correlation 
functions in Euclidean time based on the Maximum Entropy Method without a 
priori assumptions on the spectral shape. The thermal modifications of hadron 
properties are investigated in the scalar, pseudo-scalar, vector and axial-vector 
channels near and above the deconfinement transition temperature for charmonium 
systems as well as for light quarks. Moreover, the reconstructed vector meson 
spectral function allows to extract the thermal cross section for the production of 
dilepton pairs at vanishing momentum. 



1. Introduction 

Precise lattice QCD calculations of hadron masses and decay constants at 
zero temperature were performed in the recent years. A challenge is the 
investigation of modifications of hadronic properties at finite temperature, 
in particular the changes of the spectral shape. The spectral functions are 
directly relevant for experimental annihilation cross sections in relativistic 
heavy- ion collisions, e.g. the temperature dependence of the vector meson 
mass and width is related to changes of the thermal dilepton spectrum. 

On the lattice the particle masses are usually obtained from the ex- 
ponential decrease of the correlation function at zero momentum at large 
Euclidean times r: 

G(r) = fd 3 x { J H (f,r)4(0,0) ) ~ e~ m \ (1) 

where the hadronic current Jh(t, x) — ip(r, x)Yij)(T, x) is given in the scalar, 
pseudo-scalar, vector or axial- vector channel with Y £ [1,75,7^,7^75 1. At 
finite temperature T the Euclidean time extent is restricted by 1/T ren- 
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dcring a reliable determination of the particle masses increasingly difficult 
with rising temperature. 

Another possibility to extract the masses and decay constants is pro- 
vided by the hadronic spectral spectral functions an, which are related to 
the correlation functions in Euclidean time by 

oo 

G(t,p) = J d 3 xe i r* (Jh(t,x)jI(0M= J do; K(t,lu) a H {w,p,T),(2) 

o 

where K is the integration kernel in the continuum 

coshKr - 1/2T)) 
K (T '" )= sinh(c/2T) • (3) 

The reconstruction of a smooth spectral function given only a limited num- 
ber of points of the correlator in Euclidean time can be classified as typical 
ill-posed problem. Previous analyses have been done under restricting as- 
sumptions on the spectral shape. The modern solution is provided by the 
Maximum Entropy Method (MEM), which has been applied successfully to 
many similar problems in physics 1 . This approach has been first applied in 
lattice QCD by 2 . It allows to obtain the most probable spectral function 
without a priori assumptions on the spectral shape by maximizing the prior 
probability in addition to usual least-square minimization. 



2. Meson spectral functions at T — and T — oo 

Before we proceed to hadronic spectral functions at finite temperature, 
where very little is known about the spectral shape so far, the behavior 
in the limiting cases T = and T = oo should be investigated. At zero 
temperature 8— function like peaks are expected for the ground and ex- 
cited states at low energies, while the meson spectral functions should be 
proportional to lo 2 at high energies. The precise quenched QCD data of 
the CP-PACS Collaboration has been used for a detailed MEM analysis of 
ground and excited state meson masses as well as decay constants 3 . They 
have investigated pseudo-scalar and vector meson correlators on very large 
lattices (32 3 x 56 to 64 3 x 112) to at four different values of the lattice 
spacing a. The pseudo-scalar spectral function p(u) = aps/uJ 2 in figure 
1 shows several peaks for the ground and excited states. On smaller lat- 
tices only the ground state 4 and one excited state 2 can be resolved by 
the reconstructed spectral function. In general the peaks are broader and 
reduced in height for smaller lattice sizes and lower statistics. 
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The ground state masses extracted from the spectral functions are the 
same as the ones obtained from exponential fits of the correlator, while the 
mass of the first excited state could be determined much more precisely 
than before. The broad state slightly below uja — 2, similarly observed 
in 4 ' 2 , diverges in the continuum limit (fig. 1 right) and thus has been 
identified as unphysical lattice artifact (cut-off effect) of the heavy fermion 
doublers appearing in the Wilson type fermion formulations. 
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Figure 1. Pseudo-scalar spectral function at T = (left) and continuum extrapolation 
(right) taken from 3 , results for the vector meson are similar. 



At finite temperature excited states and continuum-like contributions 
are expected to gain in influence on the spectral shape. The spectral func- 
tions at T = oo for free massless quarks are known from leading order 
perturbation theory, e.g. cr{?g e = 3/8tt~ 2 lj 2 tanh(w/4T). It is manda- 
tory to verify that MEM allows to reconstruct such a continuum-like shape 
correctly. It turned out 5 that this shape cannot be obtained using the 
continuum integration kernel of eq.(3) due to large contributions of heavy 
Wilson fermion doublers at high energies. The application of the finite 
lattice approximation of the continuum kernel 



K lat (T,uj) = 2u/N T 



N T — 1 

E 

n=0 



exp(- 



-) 



4 sin 2 (u>„/2) +uj 2 



(4) 



absorbs the cut-off effects efficiently and allows an almost perfect recon- 
struction already for 8-16 points in Euclidean time. The MEM analysis 
reported in the following sections is therefore performed using this lattice 
adapted integration kernel. 

In a very recent paper 6 the authors studied the cut-off effects in detail 
in the infinite temperature limit, where the spectral functions can be cal- 
culated analytically 7 . They found considerable cut-off effects for the usual 
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Wilson fcrmion formulation, which are strongly suppressed and shifted to 
the very high energy regime using a truncated perfect fermion action. 



3. Medium modifications for light mesons 

The thermal changes of meson spectral functions in the light quark sector 
have been investigated in the range 0.4-3.0 times the chiral phase transition 
temperature. The simulations were performed in quenched QCD with a 
non-perturbatively improved fermion action on large isotropic lattices up 
to 64 3 x 24. Calculations below T c were carried out at non-vanishing quark 
mass corresponding to a pion mass of about 500 MeV. Above the critical 
temperature the simulations could be performed directly at approximately 
zero quark mass, since massless Goldstone modes are no longer present. 





Figure 2. Light pseudo-scalar (top) and vector meson (bottom) correlation (left) and 
spectral functions (right) in the range 0.4 to 3.0 T c taken from 8 . Points of the vector 
meson correlators are slightly displaced for better visibility. Only the low energy part of 
the spectral functions is given, which is not affected by cut-off effects 9 . 



Figure 2 shows the modifications of the zero momentum correlation 
and reconstructed spectral functions in the investigated temperature range. 
While crossing the phase transition temperature the change of the pseudo- 
scalar correlator is clearly obvious. Above the critical temperature the 
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correlation function slowly approaches the free meson correlator indicated 
as solid line, while the vector meson correlator is near the free curve for 
all temperatures. The corresponding spectral functions show a pronounced 
ground state peak below T c , which is transformed into a broad bump at 
higher energies above T c . In fact, this remnant of the ground state peak 
scales with the temperature and thus its position remains almost the same 
when plotted in units oiu/T instead of u/T c . 

Moreover, a degeneracy of the pseudo-scalar and scalar iso- vector corre- 
lation and spectral functions could be observed above the critical temper- 
ature indicating effective chiral Ua(1) symmetry restoration. 

A more quantitative picture of the deviations from the free quark be- 
havior can be obtained by taking the ratio of the correlator with the corre- 
sponding free meson correlator for the same lattice size, which is illustrated 
in figure 3. In the pseudo-scalar channel one can still observe a considerable 
deviation from the propagation of a free qq pair even at a temperature as 
high as 3 T c indicating a remnant of the strong interaction between the 
quarks above the chiral phase transition. The situation is different in the 
vector channel. The correlator ratio visualizes that the vector meson corre- 
lation function shows only a 10 % deviation from the free quark behavior. 
Furthermore, the value of GV(1/2T, T)/T 3 is finite, which implies that the 
spectral function ay has to vanish in the limit uo — > 0. The determination 
of the slope of the spectral function at u> — would be of major interest, 
since finite values for the transport coefficients of the QGP would require 
a ~ to in this limit 11 . Unfortunately, simulations on even larger lattices 
would be needed to increase the sensitivity of the reconstructed spectral 
functions in the very low energy region. 
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Figure 3. Correlator ratios G(tT) /G f ree (rT) for the pseudo-scalar (left) and vector 
meson (right) taken from 10 . 
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4. Thermal dilepton rates 

The thermal dilepton spectrum is accessible in heavy ion collisions and pro- 
vides an important observable to study thermal properties of the medium 
at high temperature and density 12 ' 13 . The non-perturbative in-medium 
modifications of the qq interactions are expected to influence thermally in- 
duced changes of the dilepton spectrum at low energies 14 . The differential 
dilepton rate in two-flavor QCD is directly related to the spectral function 
in the vector channel 
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a v (u,p,T) . 



(5) 



The left hand side of figure 4 gives a detailed error estimate on the vec- 
tor meson spectral function which is used to calculate the dilepton rate. 
The dashed band illustrates solely the statistical error, while the vertical 
lines illustrate the uncertainty of the MEM reconstruction in the given uj/T 
ranges. Moreover, calculations with different lattice spacings and volumes 
have been performed to ensure that the spectral shape of oy is not in- 
fluenced by cut-off or finite volume effects in the given energy regime 9 . 
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Figure 4. Vector meson spectral function (left) and thermal dilepton rate at p = 
(right) taken from 1() . 8 . 



The thermal cross section for the production of dilepton pairs at van- 
ishing momentum (right side of fig. 4) is obtained by applying eq.(5) for 
the vector meson spectral function. The leading order perturbative result 
(Born rate) and a calculation in HTL-resummed perturbation theory 15 ' 16 
are given for comparison. A small enhancement over the Born rate in the 
range 4-8 uj/T is visible for the lattice gauge theory (LGT) result at both 
temperatures, whereas the rate decreases rapidly at smaller energies. This 
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behavior is in clear contrast to the hard thermal loop calculations, which 
predict a divergent result in the infrared limit. 

It is evident that perturbative as well as lattice calculations should in- 
vestigate the low energy regime in more detail. An obvious step in this 
direction is to try to enhance the sensitivity of the statistical analysis of 
thermal correlation functions in the low energy region. A first attempt 
along this line has recently been tested in a calculation of the electrical 
conductivity and soft photon production rate of the QCD plasma 17 . 



5. Charmonia 

For the heavy cc bound states the situation is different from the light quark 
sector. Potential model calculations predict a sequential dissolution, where 
higher excitations like Xc,o [ 3 Po] and Xc,i [ 3 Pi] dissolve earlier, while the 
s-wave states J/ip [ 3 Si] and r] c I 1 So] may survive after the deconfinement 
phase transition 18 . 

Three different lattice QCD groups have investigated charmonium spec- 
tral functions so far. The simulations have been performed on large 
isotropic 19 and anisotropic lattices 20 ' 21 covering the temperature range 
0.9 to 1.9 T c . The thermal modifications of the charmonium spectral func- 
tions found by 19 are displayed in figure 5 for the pseudo-scalar (j] c and 
similar the vector meson J/ip) and for the axial- vector channel (xc.i and 
similar for the scalar Xc,o) at 0.93 and 1.25 T c . The left figure shows that 
the pronounced ground state peaks for the J/ip and r\ c survive even at 1.25 
T c , while a drastic change can be observed in the spectral functions of the 
orbitally excited states \c,o and Xc,i- The second broad bump at higher 
energies in both pictures is most likely a cut-off effect as discussed above. 




Figure 5. Thermal changes of charmonium spectral functions in the pseudo-scalar (left) 
and axial- vector channel (right) taken from 19 . 
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The results on anisotropic lattices support this finding of sequential 
dissolution. The survival of sharp ground state peaks at 1.1 T c has been 
reported for the pseudo-scalar (t] c ) and vector channel (J/tp) 20 , while a 
near degeneracy of the spectral functions in all channels has been observed 
at 1.9 T c indicating chiral symmetry restoration . Combining the results 
of all groups one may conclude that the orbitally excited states dissolve in 
the vicinity of the deconfincmcnt phase transition temperature, while the 
J/tp and T] c persist as bound states up to 1.25 T c and probably dissolve in 
the range 1.5 - 1.9 T c . 



6. Conclusions & Outlook 

The application of the Maximum Entropy Method for the reconstruction 
of spectral functions from correlation functions in Euclidean time without 
a priori assumptions on the spectral shape initiated the recent progress of 
lattice investigations in thermal physics. Already at zero temperature this 
approach leads to more precise predictions of the masses of excited states, 
but the most important point is that the study of hadronic spectral func- 
tions at finite temperature became feasible. New insight could be obtained 
in the in-medium modifications of hadronic properties in the vicinity and 
above the phase transition temperature. 

Deviations from the free quark-antiquark propagation could be observed 
in the pseudo-scalar channel up to 3 T c , while the vector meson correlation 
and spectral function are much closer to the free quark behavior. Above 
the critical temperature an effective J/a(1) symmetry restoration between 
the pseudo-scalar and scalar iso-vector states has been found. The thermal 
dilcpton production rate could be calculated for the first time from lattice 
QCD data using the spectral function in the vector meson channel. Finally, 
the pattern of sequential dissolution of the heavy charmonium states is 
supported by the spectral analysis. 

Further investigations will be certainly needed in the vicinity of the 
critical temperature and in the low energy region. All the lattice QCD 
calculations reported in this summary were performed in the quenched ap- 
proximation, which neglects the influence of virtual quark loops. The ex- 
tension to simulations with dynamical quarks will be the major challenge 
in the future. 
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